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We study the relativistic orbit of binary black holes in systems with small mass ratio. The trajectory 
of the smaller object (another black hole or a neutron star), represented as a particle, is determined 
by the geodesic equation on the perturbed massive black hole spacetime. The particle itself generates 
the gravitational perturbations leading to a problem that needs regularization. Here we study 
perturbations around a Schwarzschild black hole using Moncrief's gauge invariant formalism. We 
decompose the perturbations into I— multipoles to show that all £— metric coefficients are C° at the 
location of the particle. Summing over I, to reconstruct the full metric, gives a formally divergent 
result. We succeed in bringing this sum to a generalized Riemann's £— function regularization scheme 
and show that this is tantamount to subtract the i — > 00 piece to each multipole. We explicitly carry 
out this regularization and numerically compute the first order geodesies. Application of this method 
to general orbits around rotating black holes would generate accurate templates for gravitational 
wave laser interferometric detectors. 
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I. INTRODUCTION 

The computation of the gravitational radiation gen- 
erated by binary black holes is of great theoretical and 
observational interest. On one hand these binary systems 
are among the best candidates to be detected by the new 
generation of gravitational wave detectors. On the other 
hand the theoretically very interesting two body prob- 
lem in General Relativity remains an unsolved task for 
full numerical approaches. An appropriate astrophysical 
model to compute gravitational radiation coming from 
the capture of stars by massive black holes at the cen- 
ter of galaxies is provided by the perturbative approach 
applied to binary systems where one black hole is much 
more massive than the other. In the perturbative regime, 
linearized Einstein equations can be brought to two sim- 
ple wave equations for the two polarizations of the gravi- 
tational field [jl]j^] . This problem has been recently revis- 
ited to include initial data into the formulation in order 
to be able to start numerical integrations from finite sep- 
arations of the holes || . It also proved to be an excellent 
arena where to test the accuracy of the Longitudinal- 
conformally flat ansatz to solve the initial value problem 
for binary black holes HQ . In order to answer questions 
like what is the the displacement of the innermost sta- 
ble circular orbit (ISCO) and the rate of increase of the 
characteristic gravitational frequency due to inspiralling 
orbits, it is important to go beyond the leading approx- 
imation. Technically, one has to compute the geodesic 
trajectory of a particle in the perturbed black hole space- 
time generated by the particle itself. Taking this effect 



into account will make otherwise particle's bounded or- 
bits actually to inspiral towards the bigger black hole. 
This consistent description of the first order perturbative 
approach opens the door to second order perturbative 
studies |g] which allow to compute gravitational radia- 
tion from binary black holes with much higher accuracy 
and for systems with not so small mass ratio. Such stud- 
ies would teach us about new nonlinear physical effects 
as well as to produce accurate templates to analyze the 
forthcoming data from ground (and space) based laser 
interferometric detectors. 

Recent work on gravitational radiation reaction |||7]] 
prompted a renewed theoretical interest on this problem. 
The approach we will develop here privileges its direct 
computational implementability. Throughout this paper 
we will use results and techniques of Ref. There it 
was computed the gravitational radiation generated by 
the collision of two non-spinning black holes, one much 
less massive than the other, starting from rest at a finite 
distance. The extreme mass ratio allowed to describe the 
problem of gravitational radiation as perturbations about 
the Schwarzschild metric. Hence, instead of working with 
all ten metric perturbations, seven even parity and three 
odd parity (that identically vanish for our axially sym- 
metric problem) , the relevant perturbative information is 
organized into (one in the even parity case, two in gen- 
eral) the Moncrief waveform || 
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where we have used Zerilli's [El normalization for i\n and 



1 



notation for A = (£ + 2)(£- l)/2 . 

One of the advantages of working with ip£ is its gauge 
invariance under first order diffeomorfisms. This al- 
lows us to choose any convenient gauge, like the Regge- 
Wheeler gauge pj], to make computations, ipi satisfies 
a single (in general two, with the odd parity case) wave 
equation 



where 



dt 2 dr* 2 
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where r* = r + 2Mln(r/2M - 1), V e is the Zerilli po- 
tential, and Si(r, t), is the source term generated by the 
small hole, given in Ref. |3| . The smaller hole is described 
as a point particle of proper mass mo, its stress energy 
tensor given by 
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where is the particle 4-velocity. The two dimensional 
delta function 8 2 [f2] gives the angular location of the par- 
ticle (9 p ,(j)p). Since T MI/ , and hence Sg(r,t) are already 
proportional to mo, to first perturbative order, the radial 
trajectory r p (t), follows from the geodesic equation in the 
background geometry (Schwarzschild's here). There are 
situations when one wants to know the trajectory of the 
small hole to the next order. In the computation of the 
gravitational radiation to second perturbative order one 
needs to know the source term to second order, hence 
the trajectory of the hole on the first order metric (back- 
ground plus first order perturbations). Even within first 
order perturbation theory one would also like to go fur- 
ther in order to compute secular effects like the particle's 
bounded orbit decay around a much bigger black hole. 

For the sake of simplicity we will treat here the ra- 
dial infall of a particle into a nonrotating hole. This 
problem contains many of the relevant features that oc- 
cur for more general orbits. To first perturbative order, 
the trajectory of the particle is given by a geodesic on 
the first order metric (Schwarzschild plus first order per- 
turbations). This is so because the only "forces" acting 
on the particle are gravitational. The time component 
of the particle's four-momentum, P t = mQg tt (dt / dr) = 
—m E, is no longer a conserved quantity along the tra- 
jectory. We then have to deal with the radial and 
time (here the only nontrivial ones) components of the 
geodesic equation. We combine these two equations 
into a single equation of motion for r(t) by eliminat- 
ing the dependence on the geodesic parameter. Lin- 
earization of the resulting equation and subtraction of 
the zeroth order geodesic gives the deviation of the 
trajectory Ar p (t) from the zeroth order one r p (t) = 

d t z p = -(1 - 2M/z p ) \p- M {l P 2^^ Z ° -, directly in terms 
of Schwarzschild coordinates, ready for further applica- 
tions 
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The (numerical) integration of this expression gives the 
trajectory of the particle correct to first perturbative or- 
der, (jtiq/M). Since the metric perturbations generated 
by a particle seem to involve radial delta functions (and 
derivatives of them) evaluated at the location of the par- 
ticle (see Eqs. (||)-(||) below), the first problem to face 
here is to evaluate the connection coefficients at r p . We 
will next show explicitly that the metric is actually C° at 
the location of the particle. Connection coefficients have 
a finite jump and they can be computed as the average 
of its values at r p ± e with e — * 0. While this accounts for 
the terms coming from the radial Dirac delta in Eq. (|^), 
as we will see, the angular delta needs to be explicitly 
regularized by a different method. 



II. CONTINUITY OF METRIC COEFFICIENTS 

The tt component of Einstein's equations give us the 
Hamiltonian constraint. In the Regge- Wheeler gauge 
(hi = h a = G = 0) it is given by Zerilli's § Eq. (C7a). 
Only two metric coefficients (K l and H^) appear in this 
equation and none of its time derivatives. We recall here 
that metric coefficients have an explicit multipole index I 
since we have decomposed their angular dependence into 
tensor harmonics. Consequently K l and iff are only 
functions of t and r . Taking the Regge- Wheeler gauge 
as an intermediate step, the definition of ipe, and the 
Hamiltonian constraint we can express these two metric 
coefficients in terms of ipi only (and source terms) 
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Integration over r of the Hamiltonian constraint tells 
us that the leading behavior is given by 

d r ip e ~ k U°{r - 2M)5[r - r p ]/(X + l)/(Ar + 3M). (8) 

This can be used to prove that the metric coefficients at 
the location of the particle are actually C°, by taking 
up to second derivatives and using the equation above to 
cancel derivatives of the Dirac's delta. The same C° be- 
havior at r p can be proven for Hf . we now consider the tr 
and t6 (or tip) components of the Einstein equations that 
give us the momentum constraint. In the Regge- Wheeler 
gauge they are given by Zerilli's |] Eqs. (C7b) and (C7d). 
We combine them to eliminate the dependence on Hf 
, and obtain after integration over r that d r (d t ipe) ~ 
— k U° r p (r - 2M)d r 6[r - r p ]/(X + l)/(Ar + 3M). 

From Zerilli's Eq. (C7b) and the expressions for 
dtK and dtH 2 in terms of dtipe, we find the last metric 
coefficient in the Regge- Wheeler gauge 
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This equation together with Hq = H^ (valid for head- 
on collisions), (0), and (||) give us all metric perturbations 
on the chosen hypersurface in terms only of ipt and dfipi 
(and the source). We recall here that since our case has 
axial symmetry only even waves are generated. 



III. REGULARIZATION 

In Fig |l| we plot the results of computing by means of 
Eq. (g) the metric coefficient H 2 =(l-2M/r p )y/2t + 1H% 
along the trajectory of the particle r p (t). It clearly shows 
that it is finite at each point of the trajectory. The other 
notable feature is that for large £ the curves quickly ac- 
cumulate over the £ — -> oo asymptotic curve. Since one 
has then to sum over all £ contributions, this sum clearly 
diverges. One way of regularizing this sum is to subtract 
to each mode precisely the £ — > oo contribution, and then 
verify the convergence of the remanent series. The result 
of such regularization is shown in the lower half of Fig 
|l|. The same qualitative results are found for the other 
metric coefficients H[{r p ), K £ (r p ), and Hq{t p ) = H^{r p ). 
Non-radiative multipoles £ = 0, 1 can be found analyti- 
cally and represent the mass and linear momentum con- 
tributions of the particle to the system || . The regular- 
ization just described is ambiguous up to a finite piece. 
To give a definite prescription we note that it can be 
brought into the generalized Riemann's £— function pX|] 
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FIG. 1. The upper part of the plot shows the metric coef- 
ficient H2 along the trajectory of the particle in units of mo- 
Curves with increasing I quickly tend to superpose to that for 
L — > 00. The lower part of the plot represents the renormal- 



ized metric Ho 



tio — tin 



The inset figure gives the 



^ leading decay power, a, of H 2 ren at a fiducial r p 



procedure as follows. The numerical behavior of all met- 
ric coefficients shows that they can be decomposed into 
two pieces: One that generates the finite behavior for 
£ — > 00 and the other strongly decaying for large £ (la- 
beled with a ren below). Thus, for instance, for the per- 
turbed metric component g rr we can write 
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where in the second equality we made use of the ax- 
ial symmetry of the problem and we have chosen the 
(3— parametrization motivated by the D— dimensional ex- 
tension of the conformally flat initial value problem, 
where j3 = 4 — D. When we evaluate this metric co- 
efficient at the location of the particle we find 



H 2 (t,r p (t),6 p = 0) = 

2-/5+1/2^2 C(/ J_ 1/2) l/2) + £ 
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where the Riemann's ^-function is C( a i^) — ^2tLo(^ 
b)~ a . Numerically, we observe that (3 = 1/2 in order to 
lead to the finite £ — > 00 behavior. Since the analyti- 
cally continued C-nmction gives £(0,1/2) = we must 
subtract to each multipole just the £ — > 00 piece. The 
renormalized metric coefficients that enter into the head- 
on geodesic equation (Hq ren = H^ ren ,H( ren ) vanish 
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stages of the trajectory. This is expected for a parti- 
cle starting at rest, since the gravitational radiation is 
dominated by the bremsstrahlung mechanism. Then, ra- 
diation reaction effects become more important as the 
particle approaches the maximum of the Zerilli's poten- 
tial (around r max ss 3.1M). They tend to decelerate the 
particle with respect to the zeroth order (Schwarzschild) 
geodesies. These is what one would qualitatively expect a 
priori since the system is loosing energy and momentum 
in the form of gravitational radiation. Another impor- 
tant feature is that the radiation reaction effect is low 
^—dominated (only after renormalization) , but still the 
sum over large £ give an important contribution, and one 
has to consider higher multipoles than one usually takes 
into account when computing, for instance, the total en- 
ergy radiated reaching infinity. 

A related approach to the one presented here has been 
developed independently |llj] and applied to scalar radia- 
tion. It considers the large I dependence of the tail term 
of the reaction force , and regularizes the sum over £ 
modes by subtracting the non-convergent terms plus a 
finite part to be determined. Applied to the equivalent 
case we treated here the finite part to be subtracted van- 
ishes and both procedures coincide. Besides, this method 
was implemented to study static scalar and Electric fields 
in the Schwarzschild background , leading to the cor- 
rect known expressions of the self force. It would be 
very interesting to compare the results of our procedure 
with still other methods of regularization in order to cross 
check for possible ambiguities in the determination of the 
remaining finite parts. 



FIG. 2. The renormalized radiative piece of the reaction 
on r p . We show the regularized sum over ^ < 10 multipole 
contributions to C, as defined in Eq. (^J). Below the first order 
trajectory r\, , for mo = 0.1M, is compared to the zeroth 
order one, rp. 

for large £ as ~ £~ 3 , as we numerically roughly estimated 
(cf. Fig. [I]). This implies that the regularized connection 
coefficients scale as ~ £~ 2 for large £ and that its sum, 
to compute C, up to a finite maximum L (as is done in 
practice) scale as ~ L~ l , This can be used to estimate 
the error produced by a truncation of the series. 

From the metric we build up the connection coefficients 
that enter into the geodesic equation. The piece of the 
relative acceleration that is summed over £, denoted by C 
in Eq. (|^), needs to be regularized in the same way met- 
ric coefficients do. This piece summed over the £ < 10 
radiative multipoles is plotted in Fig. |] (t = and £ = 1 
multipoles have been solved analytically using the addi- 
tional gauge freedom to set H2 = K and, in the £ = 
case, also Hi = 0). We find the radiation reaction ef- 
fects can be qualitatively described in two regimes. First, 
the renormalized acceleration is very small in the early 



IV. DISCUSSION 

Here we reported on a first important step to im- 
prove our ability to compute gravitational radiation from 
binary black holes. The next step is to apply the 
(■—function regularization method to general bounded or- 
bits in the Schwarzschild background and compare its 
results with the energy and momentum balance esti- 
mates. Besides, our renormalization procedure is explic- 
itly gauge invariant since we use Moncrief's waveform. 
Another key problem to attack is the orbiting particle in 
a rotating black hole background. Perturbations of the 
Kerr metric are described by the Teukolsky equation and 
can be decomposed into multipoles (via spin-weighted 
spheroidal harmonics) in the frequency domain, i.e. af- 
ter a Fourier decomposition of the time dependence. This 
property would still allow to apply our regularization 
scheme and compute the decay of bounded orbits consis- 
tently to first order. Besides, it would open the possibil- 
ity to study second order perturbations Js) and obtain a 
remarkable improvement in our ability to compute gravi- 
tational radiation from binary systems with not so small 
mass ratio. This will not only be of use for LISA's detec- 
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tion of gravitational radiation from black holes in the cen- 
ter of nearby galaxies, but also relevant for ground based 
interferometers, sensible to frequencies corresponding to 
black hole / neutron star binaries. 
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